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Abstract 

Work  on  the  grant  has  concentrated  on  two  aspects:  1)  optimal 
maneuvering  of  flexible  spacecraft  and  2)  control  of  damped  distributed 
structures. 

The  problem  of  simultaneous  maneuver  and  vibration  control  of 
flexible  spacecraft  is  solved  by  means  of  a  perturbation  approach  where¬ 
by  the  slewing  of  the  spacecraft  regarded  as  rigid  represents  the  zero- 
order  problem  and  the  control  of  vibration,  as  well  as  perturbations 
from  the  rigid-body  maneuver,  represents  the  first-order  problem.  The 
perturbation  approach  to  maneuvering  of  flexible  spacecraft  was 
developed  earlier  by  the  author  of  this  report.  However,  the  present 
control  design  is  different  from  the  earlier  design,  particularly  the 
feedback  control  of  the  perturbations.  The  zero-order  problem  is  solved 
according  to  a  minimum-time  policy,  resulting  in  bang-bang  control.  The 
control  of  the  elastic  vibration  and  of  the  perturbations  from  the 
rigid-body  maneuver  is  carried  out  simultaneously  with  the  "rigid-body" 
slewing,  rather  than  sequentially,  so  that  the  overall  result  is  very 
close  to  an  ideal  minimum-time  solution.  The  maneuver  control  is  open- 
loop.  On  the  other  hand,  the  control  policy  for  the  perturbations 
involves  the  solution  of  a  time-varying  linear  regulator  problem  capable 
of  accommodating  the  disturbances  from  the  rigid-body  maneuver.  An 
exponential  factor  in  the  performance  index  forces  the  error  to  zero 
within  the  minimum-time  maneuver  period.  Feedback  control  of  the 
perturbations  is  essential  in  view  of  possible  inaccuracies  in  the 
system  parameters  affecting  the  natural  frequencies  and  mode  shapes. 

The  work  is  described  in  Ref,  1. 

Undamped  distributed  structures  represent  self-adjoint  systems, 

/ 


admitting  real  eigenvalue  and  real  orthogonal  eigenfunctions.  Control 
of  self-adjoint  systems  can  be  carried  out  conveniently  by  modal 
control.  Distributed  structures  with  proportional  damping  possess  the 
same  eigenfunctions  as  the  corresponding  undamped  structures,  so  that 
the  same  modal  approach  can  be  used  in  this  case  as  well.  Nonpropor¬ 
tional  damping  tends  to  destroy  the  self-adjointness  of  the  system,  so 
that  modal  control  is  not  as  convenient  as  for  undamped  structures.  If 
damping  is  relatively  small,  however,  it  is  possible  to  base  the  control 
design  on  the  self-adjoint  system  and  still  obtain  satisfactory  control 
performance  (Ref.  2). 
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Optimal  Vibration  Control  of  a  Flexible  Spacecraft 
During  a  Minimum-Time  Maneuver 

L.  Meirovltch  &  Y.  Sharony 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 


Abstract 


This  paper  Is  concerned  with  the  simultaneous  maneuver  and  vibra¬ 
tion  control  of  a  flexible  spacecraft.  The  problem  Is  solved  by  means 
of  a  perturbation  approach  whereby  the  slewing  of  the  spacecraft  re¬ 
garded  as  rigid  represents  the  zero-order  problem  and  the  control  of 
vibration,  as  well  as  of  perturbations  from  the  rigid-body  maneuver, 
represents  the  first-order  problem.  The  zero-order  control  is  to  be 
carried  out  in  minimum  time,  which  Implies  bang-bang  control.  On  the 
other  hand,  the  first-order  control  is  a  time-dependent  linear  quadratic 
regulator  including  integral  feedback  and  prescribed  convergence  rate. 


1.  Introduction 


The  problem  of  slewing  a  flexible  spacecraft  is  described  by  a 
hybrid  set  of  equations,  in  the  sense  that  the  rigid-body  motions  of  the 
spacecraft  are  described  by  ordinary  differential  equations,  generally 
nonlinear,  and  the  elastic  motion  of  the  flexible  parts  by  partial  dif¬ 
ferential  equations.  Practical  considerations  dictate  that  the  partial 
differential  equations  be  discretized  In  space,  resulting  in  a  set  of 
nonlinear  ordinary  differential  equations  of  relatively  high  order. 

The  problem  has  been  considered  by  Qunkins  and  Turner  (Ref.  I), 
Turner  and  Chun  (Ref.  2)  and  Chun,  Turner  and  Juang  (Ref.  3).  The 
approach  of  Refs.  1-3  was  to  minimize  a  quadratic  performance  index  for 
the  nonlinear  discretized  system.  The  problem  was  characterized  by 
prescribed  end  position  and  fixed  final  time.  A  different  approach, 
suggested  by  Ben-Asher,  Burns  and  Cliff  (Ref.  4)  and  by  Thompson, 

Junkins  and  Vadall  (Ref.  5),  consists  of  letting  the  final  time  be  free 
and  solving  a  minimum-time  problem  for  the  nonlinear  model.  Both 
approaches  involve  the  solution  of  a  nonlinear  two-point  boundary-value 
problem  for  a  high-order  system.  In  the  second  approach,  the  control  is 
open- loop. 

A  different  approach  was  developed  by  Meirovltch  and  Quinn  (Refs.  6 
and  7).  The  approach  is  based  on  the  perturbation  concept  for  solving 
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nonlinear  differential  equations.  The  zero-order  problem  consists  of 
the  rigid-body  maneuver  of  the  flexible  spacecraft,  and  is  described  by 
six  nonlinear  ordinary  differential  equations.  The  first-order  problem 
is  obtained  by  linearizing  the  system  of  equations  about  the  trajectory 
describing  the  maneuver.  The  perturbation  equations  represent  a  high- 
order  set  of  linear  time-varying  ordinary  differential  equations.  The 
rigid-body  maneuver  is  open-loop  and  the  control  of  the  perturbed  system 
is  closed-loop.  Whereas  the  feedback  control  is  carried  out  during  the 
maneuver,  the  bulk  of  the  vibration  suppression  takes  place  after  the 
termination  of  the  maneuver. 

This  paper  adopts  the  perturbation  approach  of  Refs.  6  and  7,  but 
the  control  design  is  different.,  particularly  the  feedback  control  of 
the  perturbations.  The  zero-order  problem  is  solved  according  to  a 
minimum-time  policy,  resulting  in  bang-bang  control.  The  control  of  the 
elastic  vibration  and  of  the  perturbations  from  the  rigid-body  maneuver 
is  carried  out  simultaneously  with  the  "rigid-body"  slewing,  so  that  the 
overall  result  is  very  close  to  an  ideal  minimum-time  motion.  The 
maneuver  control  is  open-loop.  On  the  other  hand,  the  control  policy 
for  the  perturbations  involves  the  solution  of  a  time-varlng  linear 
regulator  problem  (LQR)  capable  of  accommodating  the  disturbances  caused 
by  the  rigid-body  maneuver.  An  exponential  factor  in  the  performance 
index  forces  the  error  to  zero  within  the  minimum-time  maneuver  period. 
Feedback  control  for  the  perturbations  is  essential  In  view  of  possible 
inaccuracies  in  the  system  parameters  affecting  the  natural  frequencies 
and  mode  shapes. 


2.  The  Equations  of  Motion 

We  consider  a  flexible  spacecraft  consisting  of  a  rigid  hub  and  a 
flexible  appendage,  as  shown  in  Fig.  1.  For  convenience,  we  introduce 
an  inertial  reference  frame  XYZ  and  a  set  of  axes  xyz  embedded  in  the 
spacecraft  in  undeformed  state,  so  that  x  coincides  with  the  axis  of  the 
undeformed  appendage.  We  shall  refer  to  xyz  as  body  axes.  The  origin  0 
of  the  body  axes  coincides  with  mass  center  of  the  spacecraft  in  unde¬ 
formed  state.  The  general  motion  of  the  spacecraft  can  be  described  in 
terms  of  the  translation  and  rotation  of  the  body  axes  relative  to  the 
inertial  frame  and  the  elastic  motion  of  the  appendage  relative  to  the 
body  axes.  To  this  end,  we  denote  by  R  the  position  vector  of  0  rela¬ 
tive  to  the  inertial  space,  by  r  the  nominal  position  vector  of  a  point 
in  the  spacecraft  relative  to  the  body  axes  and  by  u(r,t)  the  elastic 

displacement  vector  of  a  typical  point  in  the  appendage.  Moreover,  we 
denote  by  «(t)  the  angular  velocity  of  the  body  axes  relative  to  the 
inertial  space.  The  equations  of  motion  represent  a  hybrid  set  of  dif¬ 
ferential  equations  consisting  of  six  ordinary  differential  equations 
for  the  rigid-body  motion  of  the  body  axes  and  three  partial  differen¬ 
tial  equations  for  the  elastic  motion  relative  to  the  body  axes.  The 
partial  differential  equations  for  u  can  be  replaced  by  3N  ordinary  dif¬ 
ferential  equations,  where  N  is  the”number  of  degrees  of  freedom  used  to 
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represent  each  component  of  u.  Reference  6  presents  the  complete  set  of 

equations.  Our  Interest  Is  In  a  special  type  of  motion,  so  that  the 
general  equations  are  not  listed  here. 

Let  us  consider  the  single-axis  maneuver  of  the  spacecraft,  and  In 
particular  the  slewing  of  the  spacecraft  In  the  yz-plane,  which  amounts 
to  the  rotation  about  the  x-axIs  through  a  given  angle.  Ideally,  the 
maneuver  Is  to  be  carried  out  as  If  the  spacecraft  were  rigid.  In 
practice,  a  pure  rigid-body  maneuver  Is  difficult  to  achieve,  so  that 
the  maneuver  will  excite  elastic  motions,  which  In  turn  will  cause  the 
spacecraft  to  deviate  from  the  rigid-body  maneuver.  Due  to  the  con¬ 
figuration  of  the  spacecraft  under  consideration  (see  Fig.  1),  the 
equations  for  the  rigid-body  and  elastic  displacements  in  the  yz-plane 
and  the  rotation  about  the  x-axIs  are  decoupled  from  the  remaining 
equations  of  motion.  Because  these  equations  correspond  to  the  very 
same  motions  involved  In  the  slewing  maneuver,  we  concentrate  on  these 
equations  alone.  Introducing  the  necessary  simplifications  In  Eqs.  (18) 
of  Ref.  6,  the  equations  are 

mR  +  c2$Tg  +  2a  Cj£Tg  -  (fl^  +  tt*c2)$Tg  =  CTF  (la) 

Ic«x  +  IJT9  ♦  cos  «  Ry$Tg  =  M  (lb) 


ma9  +  V  +  (KA  *  nzMA)9  +  ^2?  *  9 

where  m  is  the  total  mass  of  the  spacecraft,  R  =  (Ry  RZ)T  is  the 
position  vector  of  the  mass  center  and 


(lc) 


(2) 


is  the  rotation  matrix  from  axes  XYZ  to  axes  xyz,  in  which  a  is  the 
maneuver  angle.  Of  course,  nx  =  a,  a*  =  a  are  the  maneuver  angular 
velocity  and  acceleration,  respectively.  In  addition,  u(r,t)  = 

[0  uz(y,t)]T  Is  the  elastic  displacement  vector,  where 

u2(y.t)  »  »T(y)g(t)  (3) 

in  which  $(y)  is  a  vector  of  admissible  functions  and  q(t)  Is  a  vector 
of  generalized  coordinates.  Moreover, 

♦  -  *  dmA,  i>  =  Jm  y$  dmA  (4a, b) 

A  A 

where  mA  is  the  mass  of  the  appendage.  Other  quantities  entering  into 
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Eqs.  (1)  are  the  total  mass  moment  of  Inertia  Ic  of  the  undeformed 
spacecraft  about  x  and  the  appendage  mass  and  stiffness  matrices 

ma  =  J*  $(y)$t(y)  dmA*  ka  3  t***1  (5a,b) 

A 

where  (  ,  ]  represents  an  energy  inner  product  (Ref.  8). 

Finally,  we  assume  that  there  are  two  force  actuators  Fy  and  Fz  and 

one  torquer  Mx  acting  at  the  mass  center  of  the  rigid  body  and  p  torque 

actuators  MAj  acting  on  the  elastic  appendage  at  the  points  y  =  y^  (1  = 

1,2,..., p).  The  torque  actuators  acting  on  the  appendage  can  be  expres¬ 
sed  as  a  distributed  actuator  torque  in  the  form 

P 

%  ■  ^  -  »1>  (6) 

where  «(y  -  y^)  are  spatial  Oirac  delta  functions.  Then,  the  forcing 
terms  appearing  In  Eqs.  (1)  can  be  written  as  follows: 

P 

F-  m  Fyl  +  F£-  M  '  Mx  +  A>a  "Ax  dDA  *  Mx  +  ^  MA1  (7s.») 

P 

9  =  /da  raAxi'  d0A*1I1MA1«’(7f)*E?!A  <70 

where  E  =  ( $ ' (y^)  ♦ * (y2 )  •••  $'(yp)l  an  NxP  modal  participation 

matrix,  in  which  primes  denote  derivatives  with  respect  to  y,  and  M.  = 

T 

[MAi  MA2  *•*  MAp^  a  P~vector  actuator  torques  acting  on  the 

elastic  appendage.  The  reason  for  using  torquers  Instead  of  thrusters 
on  the  elastic  appendgee  will  become  obvious  later. 


The  Perturbation  Approach  to  the  Maneuvering  Problem 


We  consider  the  case  in  which  the  maneuvering  dynamics  can  be 
assumed  to  consist  of  some  large  terms  associated  with  the  ideal  rigid- 
body  maneuvering  and  some  small  terms  associated  with  the  elastic 
motions  and  the  perturbations  they  cause  In  the  rigid-body  motions. 
Consistent  with  terminology  used  In  perturbation  analysis,  we  refer  to 
the  large  terms  as  zero-order  terms  and  to  the  small  terms  as  first- 
order  terms,  and  denote  them  by  subscripts  0  and  1,  respectively. 

Hence,  we  write 

?  =  ?0  +  C0-P  a  3  a0  +  °1  (8atb) 
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where  Rj  Is  expressed  In  terms  of  components  along  the  body  axes,  Rj  = 
(Ry  Rzj.  Note  that  CQ  =  C(aQ).  Moreover,  g  is  regarded  as  a  first- 
order  quantity.  It  can  be  verified  that 


So 


♦  cI(R 


0V~1  +  ao?l^ 


(9a) 


B  *  ?o  +  co(Bi  +  ao~i  +  2aoBi  +  ao~i^ 

+  (aQP  -  OqORj) 


Bq  +  colBi  +  2aop?i 


(9b) 


a 


x 


where  aQ 


fl0  +  “1*  hx  =  60  +  al 
OqP,  in  which 


P  = 


0 

1 


-1 

0 


(9c,d) 


(10) 


is  a  permutation  matrix.  Similarly,  the  forcing  terms  can  be  divided  as 
follows: 


P0 


Pl’ 


M  =  Mrt  +  M 


1 


whereas  Q  is  assumed  to  be  of  first  order. 


(11a, b) 


Inserting  Eqs.  (8)-( 11)  into  Eqs.  (1)  and  neglecting  second-order 
terms  in  the  perturbations,  we  obtain  the  zero-order  equations 


mR0  -  CqFq  » 


^cC0  =  M0 


(12a, b) 


where,  considering  Eqs.  (7),  FQ  =  FyQj  +  F^k,  MQ  *  Mx(J.  From  Eqs. 

(12),  it  Is  obvious  that,  owing  to  the  fact  that  the  motion  is  referred 
to  the  mass  center  C,  the  translational  and  rotational  equations  are 
independent  of  each  other.  Because  the  Interest  lies  in  a  slewing 
maneuver,  there  is  no  loss  of  generality  in  assuming  that  RQ  = 

Rq  *  Rg  *  0.  This  permits  us  to  dispense  with  Eq.  (12a),  so  that  we  are 

left  with  a  single  zero-order  equation,  Eq.  (12b).  Moreover,  we  obtain 
the  first-order  equations,  which  can  be  written  in  the  compact  form 


Mx  +  Gx  +  Kx  =  X 


where 


x  =  [R{  Bl  gT]T,  X  = 


P1 

M, 


9  -  a0* 


(13) 


(14a, b) 
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are  (3+N)-d1mens1onal  displacement  and  generalized  force  vectors, 
respectively.  In  which 


El  '  Fyli  +  f2l^-  M1  '  Mxl  *  ^ 

and  Q  Is  given  by  Eg.  (7c),  and 


A1 


ml  0  e2$T 

2mflgP  0 

•2oo?i»T 

M  = 

0  I  *T 

,'  G  = 

0  0 

0 

c 

-  T  u 

2Pq-~1  ^ 

0 

t?2  -  MA 

- 

m(Q0P  -  flgl) 

0  -(“oSl  +  q??2^T 

K  = 

0 

0 

0 

5(«oSl  •  no?2> 

0 

KA  “  n0MA 

(15a, b) 


(16a, b) 


(16c) 


are  coefficient  matrices,  where  e^  *  [1  0]  and  e2  *  [0  1]  are  stand¬ 
ard  unit  vectors.  Note  that  G  and  K  depend  on  time  through  qq  and  Qq. 

It  will  prove  convenient  to  express  the  generalized  force  vector  in 
terms  of  actual  actuator  forces  and  torques.  To  this  end,  we  Introduce 
the  notation 

Fyl  =  Ti»  Fzi  =  T2*  Mxl  =  T3*  MAi  =  T3+i ’  1  =  (17a) 


0  * 


0 

0  ... 

0  " 

i  ! 

0 

0  ... 

0 

1  ... 

1_ 

,  f  = 

°  « 

E 

J  T.T 


and  consider  Eqs.  (7c)  and  (15)  to  obtain 
X  =  0T  -  Ogf 


(17b, c) 


(18) 


where  T  is  the  (3+p) -vector  of  actuator  forces.  Moreover,  we  wish  to 
express  the  matrices  G  and  K  In  the  form 
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G  =  2QqG 


where 


Kc  +  °0G  “ 


nOKs 


0 

0 

0 


(19a. b) 


(20a) 


and  we  recognize  that  G*  Is  skew  symmetric,  K  is 
a  symmetric  matrix,  and  Ks  is  symmetric  as  well. 


(20b, c) 


the  constant  part  of  K, 


4.  Pseudo-Modal  Equations  of  Motion 


Equation  (13)  represents  a  set  of  3+N  linear  ordinary  differential 
equations  with  time-dependent  coefficients.  For  future  reference,  we 
wish  to  rewrite  the  equations  In  a  different  form.  To  this  end,  we 
observe  that  by  letting  *  nQ  *  0  Eq.  (13)  reduces  to  a  time-invariant 

set  of  equations.  The  eigenvalue  problem  for  the  time-invariant  system 


has  the  form  K  U  =  MUa,  where  U  Is  a  (3+N)x(3+N)  matrix  of  eigenvectors 
c  2  2  2 

and  a  =  dlagfO  0  0  ug  ...  Is  a  diagonal  matrix  of  eigen¬ 
values,  in  which  u^(1  *  1,2, ...,N)  are  recognized  as  the  natural  fre¬ 
quencies  of  the  nonmaneuvering  spacecraft.  Because  M  and  Kc  are  real 
and  symmetric  and,  moreover,  M  is  positive  definite  and  Kc  is  positive 

semldeflnite,  the  eigenvectors  are  real  and  orthogonal  with  respect  to  M 
and  Kc  and  the  nonzero  eigenvalues  are  real  and  positive  (Ref.  2).  The 

T  T 

eigenvectors  can  be  normalized  so  as  to  satisfy  U  MU  *  I,  U  K  U  =  a.  It 
can  also  be  verified  that 


T  * 

u'g  U 


5, 


UTKU  =  K 
s 


where 


G  = 


In  which 


K 


i  o  o  ; 

o  1  o»o 

9-9_94- 


(21a, b) 


(22a, b) 
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rN  +  UW^2Z 


and  we  note  *-;.at  IN  is  the  N*N  identity  matrix  and  U 
right  submatrix  of  U. 


T 

22 


(23) 

the  N*N  lower 


Next,  let  us  introduce  the  linear  transformation 
x(t)  *  Uv(t)  (24) 

into  Eq.  (13),  multiply  on  the  left  by  U  ^  and  obtain 

v  +  2QqGv  +  (a  +  nQG  -  «q^)v  =  V  (25) 

where 

V  =  UTX  =  UTDT  -  nQUTf  (26) 

Equation  (25)  represents  a  set  of  pseudo-modal  equations. 

Due  to  the  nature  of  the  coefficient  matrices,  G,  a  and  K,  as  well 
of  the  excitation  vector  V,  the  pseduo-modal  equations  can  be  separated 
into 

-T  +  ^n0P^T  +  («0P  ”  fl0^-T  =  ~T  =  ^1  ^2^  (27a) 

VR  =  VR  =  1.  Tj  (27b) 

J=1  J 

-E  +  ^aE  ”  n0^22^-E  =  -E  =  ^22^-E  ~  ^0^22-  (27c) 

where  =  (v^  is  a  vector  of  perturbations  In  the  rigid-body 

translations,  v^  =  V3  is  the  perturbation  in  the  rigid  body  rotation, 

-E  *  Iv4  v5  •••  v3+n1T  a  vector  of  pseudo-modal  coordinates  corres¬ 
ponding  to  elastic  motions,  a^  =  diag  (u^  •••  wjj)  and  = 

[T4  Tg  ...  T 3+p ]  is  a  vector  of  torquers  on  the  elastic  appendage. 

Equations  (27)  Indicate  that  the  rigid-body  translations,  rigid-body 
rotations  and  elastic  motions  can  be  treated  Independently.  This 
justifies  the  choice  of  actuators  in  Sec.  2. 


5.  General  Control  Policy 

The  control  can  be  divided  into  two  tasks  to  be  carried  out 
slmultaneoulsy.  The  first  task  consists  of  the  rigid-body  slewing  of 
the  structure  according  to  Eq.  (12b).  The  second  task  consists  of 
designing  a  regulator  to  suppress  the  perturbations  from  the  rigid-body 
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slewing.  To  this  end,  we  use  Eqs.  (27).  Of  course,  the  actuator  forces 
and  torques  represent  the  controls  for  both  tasks. 


1 .  Rigid-body  slewing. 

For  the  rigid-body  slewing,  we  use  a  minimum-time  control  policy, 
which  Implies  bang-bang  control  (Ref.  9).  Assuming  an  ideal  actuator 
and  recognizing  that  nQ  =  aQ,  where  aQ  is  the  rigid-body  slewing  angle, 

Eq.  (12b)  represents  a  linear  time-invariant  system  with  two  real  eigen¬ 
values.  As  a  result,  there  is  only  one  switching  time,  namely  halfway 
through  the  maneuvering  time  interval.  Hence,  denoting  by  t0  the  ini¬ 
tial  and  by  t^  the  final  time,  the  switching  time  is  simply  t^  =  (t^  - 
tg)/2,  so  that  the  control  law  is 


c  for  tQ  <  t  <t{ 
n0  -  c  for  t^  <  t  <  tf 


(28) 


Of  course,  the  torque  is  Mg  =  I^q,  so  that  c  must  be  such  that  the 

actuator  operates  at  the  saturation  point.  From  Eq.  (28),  we  conclude 
that  the  slewing  angular  velocity  is 

fct  for  tg  <  t  s  t^ 

-  c(t  -  2tj)  for  tj  <  t  <  tf  (29^ 

and  the  slewing  angle  is 


“0  =< 


^  ct2  for  tg  <  t  <  tf 
-  ct2  -  c(i  t2  -  2txt)  for  t,  <  t  <tf 


’1 


(30) 


and  we  note  that  Eqs.  (29)  and  (30)  imply  that  tg  =  0. 
ii .  Perturbation  suppression 


The  perturbations  are  governed  by  Eqs.  (27)  and  they  represent  a 
tlme-var-ving  system  subjected  to  persistent  disturbances  generated  by 
the  inertial  forces  resulting  from  the  rigid-body  slewing.  These  per¬ 
turbations  are  to  be  suppressed  during  the  maneuver,  i.e.,  during  the 
time  time  interval  tf  -  tg.  The  time-dependent  functions  in  Eqs.  (27) 

are  Ug(t)  and  Og(t),  which  represent  the  commanded  angular  velocity  and 

acceleration  of  the  "rigid  body"  slewing.  Hence,  they  are  both  known  a 
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Another  aspect  of  the  problem  Is  the  high  dimensionality.  Due  to 
the  on-line  computer  limitations,  as  well  as  Inaccuracies  In  the  model 
higher  states,  the  originally  truncated  model,  Eq.  (27c),  must  be 
truncated  again.  The  newly  truncated  model,  referred  to  as  the  reduced- 
order  model  (ROM),  will  be  controlled  by  an  optimal  reduced-order  con¬ 
troller  (ROC),  which  will  accommodate  disturbances  during  maneuver. 

The  truncation  mentioned  above  is  an  open-loop  truncation  in  which 
the  higher  states  are  ignored  on  the  assumption  that  their  excitation 
during  maneuver  is  minimal.  Moreover,  any  damping  inherent  in  the 
system  tends  to  damp  out  these  states  the  fastest.  The  validity  of  this 
assumption  can  be  verified  by  simulations. 

The  drawbacks  of  controlling  a  distributed-parameter  system  by  a 
ROC  are  the  well-known  control  and  observation  spillover  (Refs.  10  and 
11).  We  assume  in  this  paper  that  observation  spillover  effects  are 
mitigated  by  using  a  sufficient  number  of  sensors,  i.e.,  2  translational 
displacement,  1+N  angular  displacement,  2  translational  velocity  and  1+N 
angular  velocity  sensors,  where  the  sensors  are  placed  so  as  to  permit 
an  accurate  estimate  of  the  controlled^ state  vector. 

As  pointed  out  in  Sec.  4,  the  rigid-body  translations,  the  rigid- 
body  rotation  and  the  elastic  coordinates  in  the  pseudo-modal  equations, 
Eqs.  (27),  can  be  treated  independently.  In  particular,  we  first  design 
controls  for  the  elastic  perturbations.  Then,  we  design  the  control  for 
the  rotational  perturbation,  taking  into  account  the  fact  that  the  actu¬ 
ators  for  the  elastic  motion  perturb  the  rotational  motion,  as  can  be 
concluded  from  Eq.  (27b).  The  rigid-body  translations  are  entirely 
decoupled  from  the  rotational  and  elastic  perturbations.  Concentrating 
on  the  elastic  motions,  we  can  introduce  the  elastic  state  vector 

z(t)  =  [vT(t)  »  vT(t)]T  =  (zj(t)  *zJ(t)jT,  where  zQ(t)  =  (vj(t)  ;  vj(t)] 

is  the  controlled  part  of  the  state  vector  and  z^t)  =  (Vy(t)  T  Vy]  is 

the  uncontrolled  part.  Then,  considering  Eq.  (27c),  the  state  equations 
can  be  written  in  the  form 


588 


OPTIMAL  VIBRATION  CONTROL  DURING  A  MINIMUM-TIME  MANEUVER 


TJC 


__0__ 

n27 

r  true 


.i„9 

i 

»  0 
I 


•  Au 


Au  +  °0KU 


(32c, d) 


BC  = 


= 


»  ?C 


0 

ic 


,  R, 


-U 


(32e-h) 


where  we  introduced  the  notation  U^qE  »  Eg,  U^yE  =  Ey,  U^g*  =  ^c* 
B22U~  *  -U*  an<*  we  n0^e  *^a^  dimensions  °f  the  various  submatrices 
are  consistent  with  the  dimensions  of  Zg  and  Zy,  which  are  2Ng  and  2Ny, 
respectively.  Of  course,  dim  Zg  +  dim  Zy  =  2N. 


6.  Control  of  the  Elastic  Reduced-Order  Model  (ROM) 

The  elastic  ROM  is  obtained  from  Eq.  (31)  by  Ignoring  the  uncon¬ 
trolled  state  vector  Zy  and  neglecting  the  coupling  matrix  Acu.  Hence, 

we  propose  to  retain  for  control  the  reduced-order  model 


~C  *  AC?C  +  BC-E  " 


^0-C* 


yc  s  Mc?c 


(33a, b) 


where  the  matrices  Ag  and  Bc  and  the  vector  Rg  are  given  by  Eqs.  (32a), 
(32e)  and  (32g),  respectively,  and  tig  Is  given  by  Eq.  (28),  in  which  c 
Is  a  given  constant.  Moreover,  y^  Is  the  output  vector  and  It  repre¬ 


sents  the  quantity  to  be  controlled.  In  our  case,  the  output  vector  has 
two  components,  namely,  the  angular  deflection  and  velocity  at  various 
locations  y^*y2»**‘»yt  of  the  beam.  Hence,  considering  Eqs.  (3),  (14a), 

(24)  and  (27),  and  recalling  the  definition  of  Zg,  we  have 


FTU 


22 


FTU 


22 


(34) 


where  F  =  (<»'(y^)  ...  ♦'(yl)]  is  an  i  *  2Ng  matrix.  In  which 

primes  denote  derivatives  with  respect  to  y. 


The  object  of  the  control  Is  to  drive  yg  to  zero  by  the  end  of  the 
maneuver,  or 


*C(t) 


=  0 


(35) 
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Because  of  the  persistent  disturbance  -  Q0Rr.  for  Eq.  (35)  to  be  satis¬ 
fied  optimally,  we  must  also  have  U’L 

11m  BcJE(t)  *  cRc  (36) 

t-*t^ 

For  Eq.  (35)  to  be  satisfied,  the  range  space  of  Rg  must  be  contained  In 
the  range  space  of  Bc.  Hence,  because  the  order  of  the  ROM  is  equal  to 
the  dimension  of  Rg,  the  number  of  actuators  must  be  equal  to  one  half 
the  order  of  the  ROM,  or  p  =  Nr*  Of  course,  the  location  of  the 

T 

actuators  must  be  such  that  Eg  *  U^gE  has  rank* 

The  problem  of  disturbance  accommodation  using  an  LQR  was  discussed 
In  Refs.  12-17,  among  others.  In  view  of  the  dimension  of  Rg,  as  well 

as  our  goal  of  satisfying  Eq.  (35)  optimally,  the  approach  suggested  In 
Refs.  12  and  13  appears  suitable.  Moreover,  to  force  ^g(t)  to  zero 

within  the  time  Interval  (tg.tf),  we  include  an  exponential  term  In  the 

performance  index.  This  term  Is  equivalent  to  prescribing  a  degree  of 
stability  to  an  Infinite  final  time,  time-invariant  problem,  as  sug¬ 
gested  In  Ref.  17.  In  the  Appendix,  it  Is  shown  that  this  exponential 
term  does  provide  the  minimal  rate  of  convergence  even  in  time-varying 
systems. 


Let  us  consider  a  new  variable  vector  u(t)  defined  by 
Bcu(t)  *  BcTE(t)  +  Fc(t) 


where 


5c’ 


cRg  for  tg  <  t  <  tj 
cRg  for  t^  <  t  <  tf 


(37) 

(38) 


and  Introduce  the  new  state  vector  Xg(t)  =  (zj(t)  uT(t)lT.  Then,  the 
new  state  equations  can  be  written  in  the  form 


xc(t)  =  H(t)xc(t)  +  Gu^t),  u^t)  =  u(t) 
In  which 
H  = 


[AC  1  Bc 

•  G  * 

■  *■ 

_0_ 

o 

o 

1 _ 

I 

k  « 

(39a, b) 

(40a, b) 


We  propose  to  derive  an  optimal  control  by  minimizing  the 
performance  measure 

J  =  xJ(Tf)Hxc(Tf)  +  JTf  e^xjQjXg  +  ujs^)  dt  (41) 

where 
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*1- 


gj.o 

0  !  rj 


,  Q  =  mJqqMc 


(42a, b) 


The  matrices  Q,  R  and  S  are  real  symmetric  and  positive  definite.  More¬ 
over,  Tq  and  Tf  are  equal  to  tg  and  tj  for  the  first  stage  and  to  t^+ 

and  tf  for  the  second  stage. 

The  formulation  given  by  Eqs.  (39)  and  (41)  can  be  transformed  Into 
a  standard  LQR  formulation  by  defining  the  new  state  and  control  vectors 

xc(t)  =  eatxc(t),  ut(t)  =  eatu1(t)  (43a, b) 

In  terms  of  the  newly  defined  vectors,  the  state  equations  have  the  form 

xc(t)  =  H(t)xc(t)  +  Gu1(t),  x  c ( Tg )  =  e  0xc(TQ)  (44a,b) 

a 

where  H(t)  =  H(t)  +  al.  Similarly,  the  modified  performance  measure  Is 

(45) 


=  5g(T^)Hx^(T^)  +  ^Tg^C^l-C  * 


Inserting  Eqs.  (43)  Into  Eq.  (45),  It  Is  easy  to  verify  that  0  =  J. 

a 

The  optimal  control  law  can  be  shown  to  have  the  form 

11,(1)  >  u(t)  =  -  S'1GTP(t)!«c(t)  (46) 

where  P(t)  Is  the  solution  of  the  Riccatl  equation 

P(t)  =  -  P(t) (H(t)  +  all  -  [ HT ( t )  +  aljP(t)  +  P(t)GS_1GP(t)  -  (47) 

and  Is  subject  to  P(T^)  =  H.  We  can  partition  the  Riccatl  matrix  as 
follows: 

rp,,  ;  Pi,i 

(48) 


rpn_ 

P12 

pt 

1/12 

P22. 

where  Pjj  Is  2Ngx2Ng,  P^  Is  2Ng*Ng  and  P Is  Ng*Nc.  Then,  Eq.  (47) 

can  be  rewritten  In  terms  of  the  original  system,  Eqs.  (39)  and  (41),  in 
the  form 


P11  *  '  P11AC  '  ACP11 


-  2aPll  +  P12S“lpl2  -  Q 


P12  =  ‘  PUBC  ‘  ACP12  *  2aP12  +  P12S  lp22 
P22  3  ‘  P12BC  *  8CP12  “  2aP22  +  P22S  lf>22  "  R 


(49a) 

(49b) 

(49c) 
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where  the  submatrices  of  P  are  subject  to  P^fi^)  *  H^O.j  »  1,2). 

Next,  we  wish  to  express  the  actuator  vector  T^(t)  In  terms  of  the 

pseudo-modal  coordinates.  Inserting  Eqs.  (41b)  and  (48)  Into  Eq.  (46), 
we  obtain 

u(t)  *  -  S_1pJ2{t)*c(t)  -  S_1Pz2(t)u(t)  (50) 

But,  u(t)  Is  related  to  the  actuator  vector  via  Eq.  (37).  Hence,  using 
Eqs.  (33a)  and  (37),  recalling  Eqs.  (28)  and  recognizing  that  Fc(t)  *  0 
over  each  control  stage,  we  conclude  that  the  actuator  vector  satisfies 

TE(t)  *  -  S-VI2(t)  -  P22(t)B+Ac(t)|2c  -  S-1P22(t)B^ic(t)  (51) 

where  b£  »  (BqB^)  *bJ  *  [0;  (eJee)  *eJ)  *  101  E^ 1 ,  In  which  the  dagger 
denotes  the  pseudo-inverse  of  a  matrix.  Introducing  the  notation 

S_1[P{2(t)  -  P22(t)BjAc(t)l  *  (K|(t)  |  K2(t)l  (52a) 

s"lp22(t)Bc *  s"lp22(t)[0i  EC]  =  10 i  h(t)]  (52b) 

the  actuator  vector  can  be  written  in  the  form 

IE(t)  -  Te(T0)  -  IK1(»)vc(t)  .  K2(t)yc{,)  ♦  K3(,)yc(,)|dt  (53) 

where  Tq  is  either  tQ  or  tj,  depending  on  the  control  stage.  For 
convenience,  we  take  (Tq)  *  0. 

It  should  be  pointed  out  that  the  above  Integral  feedback  control 
is  capable  of  compensating  for  constant  unknown  disturbances,  such  as 

might  occur  when  the  vector  *  In  the  disturbance  vector  -  nn4i  is  not 
known  exactly.  *  u“ 


7.  Response  of  the  Full  (Discretized)  Model 

We  refer  to  the  discretized  model  given  by  Eq.  (31)  as  full, 
although  It  represents  a  truncated  model  relative  to  the  distributed 
system.  The  Idea  is  that  the  states  Ignored  In  arriving  at  Eq.  (31)  do 
not  participate  anyway.  The  closed-loop  system  Is  obtained  by  intro¬ 
ducing  Eq.  (53)  Into  Eq.  (31).  Because  the  feedback  control  law  Is  In 
Integral  form,  we  will  find  It  convenient  to  introduce  the  expanded  3N- 

state  vector  given  by  w  *  [wj  ;wJ]T  *  lyjj  yj  ;  yj;  vj;  vj  ;  vJ]T.  Then, 

using  Eqs.  (31)  and  (53),  we  obtain  the  closed-loop  equation 


w(t)  =  A(t)w(t) 


(54) 


592 


OPTIMAL  VIBRATION  CONTROL  DURING  A  MINIMUM-TIME  MANEUVER 


where 


In  which 


AC  |  ACU 

Auc«  Au 


2ao“oKC  ■  ECK1  1  nO-KC"  Ac 


I 

- L. 

c  k  I 

tc  2  ■ 


(56a) 


2ao“oKcu  !  aoKcu  ! 


(56b) 


2nonoKuc 


“  EUK2  I 


(56c) 


2n0a0KU  J  n0KU 


(56d) 


where  we  regarded  as  zero. 

In  many  structures,  the  off-diagonal  entries  of  the  matrix  K^,  Eg. 

(23),  are  very  small  compared  to  one.  This  Is  particularly  true  in  the 
case  under  Investigation.  In  such  cases,  the  coupling  matrix  Acu,  Eq. 

(56b),  Is  very  small  and  can  be  ignored  without  affecting  the  results 
very  much.  Under  these  circumstances,  we  conclude  from  Eqs.  (54)  and 
(55)  that  the  uncontrolled  states  do  not  affect  the  controlled  states. 

On  the  other  hand,  because  Au(.  *  0,  and  the  terms  preventing  Auc  from 

being  zero  can  be  traced  to  the  feedback  control,  we  conclude  that  the 
system  Is  subjected  to  control  spillover.  It  is  assumed  here  that  the 
controlled  states  are  fully  observable.  Hence,  the  control  spillover 
cannot  lead  to  instability,  although  significant  performance  degradation 
can  occur  If  the  number  of  controlled  modes  is  too  small. 
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8.  Numerical  Example 

The  structure  considered  Is  as  shown  in  Fig.  1,  In  which  the  length 
of  the  beam  Is  L  =  24  ft.  Moreover,  the  first  seven  natural  frequencies 
of  vibration  in  the  yz-plane  are  0.32,  1.8,  5.0,  9.8,  16.2,  24.2  and 
33.8  (Hz).  A  damping  factor  of  0.01  is  included. 

Two  different  maneuvers  were  considered.  The  first  consisted  of  a 
180°  slew  in  7.05s  and  the  second  of  20°  slew  in  3.98s. 

The  ROM  was  assumed  to  have  four  flexible  degrees  of  freedom  and 
the  originally  discretized  system  was  assumed  to  have  seven  flexible 
degrees  of  freedom.  There  were- seven  actuators,  three  rigid-body  actua¬ 
tors  mounted  on  the  hub  and  four  actuators  mounted  on  the  beam  at  the 

points  y.j  =  1L/4  (i  =  1,2, 3, 4).  The  values  of  the  control  parameters 

were  taken  as  Q  *  R  *  I,  S  =  0.1  I,  a  *  5.0.  Control  is  terminated  at 


The  rigid-body  maneuver  is  shown  in  Fig.  2.  It  Is  an  ideal  min¬ 
imum-time  maneuver.  Figure  3  shows  the  moment  output  of  the  actuators 
on  the  beam.  Every  actuator  generates  a  near-constant  moment,  compen¬ 
sating  for  the  constant  disturbance.  The  effect  of  the  integral  feed¬ 
back  is  obvious.  Figure  4-7  show  the  angular  error  (measured  in 
degrees)  at  the  tip  of  the  beam,  at  which  point  the  deflection  Is  the 
largest.  The  error  represents  the  difference  between  the  actual  dis¬ 
placement  and  the  rigid-body  displacement,  where  the  latter  can  be  ob¬ 
tained  from  Fig.  2.  The  characteristics  of  the  error  for  both  maneuvers 
are  similar,  although  the  first  maneuver,  being  faster,  exhibits  a 
larger  error.  Figures  4  and  6  show  the  error  of  the  controlled  model 

alone.  Convergence  to  zero  (10“^  degrees)  is  achieved  during  the 
control  interval  (Tg.Tf)  by  choosing  a  proper  convergence  factor  a.  The 

"full  model"  experienced  an  error  due  to  control  spillover  into  the 
uncontrolled  states,  as  well  as  to  uncompensated  disturbances.  This 
error  is  shown  in  Figs.  5  and  7.  The  Inherent  damping  in  the  system 
causes  this  error  to  decay  after  the  termination  of  the  control. 


9.  Conclusions 


For  maneuvers  characterized  by  angular  velocities  not  much  higer 
than  the  lowest  natural  frequency  of  the  beam,  the  time-varying  terms  in 
A^(t)  tend  to  be  negligible  in  the  case  of  closed-loop  control.  In  such 

cases,  the  solution  of  the  Riccati  equation  requires  less  computational 
effort,  as  does  not  need  updating.  Note,  however,  that  the  Riccati 

matrix  still  depends  on  time. 

The  errors  result  mainly  from  control  spillover  into  the  uncon¬ 
trolled  states.  The  contributions  to  the  errors  are  primarily  from  the 
lowest  uncontrolled  states.  The  error  can  be  reduced  by  increasing  the 
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number  of  controlled  modes.  Even  for  the  model  used  here,  the  error  is 
not  very  large.  Indeed,  the  error  In  slope  at  the  tip  of  the  beam  drops 
to  less  than  0.1*  at  t  *  1.1  tf.  This  good  performance  can  be  attrib¬ 
uted  to  the  large  separation  in  the  natural  frequencies,  which  is  typi¬ 
cal  of  one-dimensional  structures  such  as  beams.  For  two-  and  three- 
dimensional  structures,  the  natural  frequencies  tend  to  be  spaced  more 
closely,  so  that  the  error  is  likely  to  increase. 
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Appendix 

The  purpose  of  this  Appendix  is  to  examine  the  convergence  rate  of 
the  controlled  state  Xj,(t).  To  this  end,  we  recall  that  the  optimal 

control  law  for  the  system  described  by  Eqs.  (44)  and  subject  to  the 
performance  Index  (45)  Is  given  by  Eg.  (46),  In  which  P(t)  Is  the 
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solution  of  the  matrix  Rlccatl  equation,  Eq.  (47).  The  matrix  P(t)  Is 
symmetric  and  positive  definite  for  te[tQ,t^). 

To  help  with  the  developments,  we  Introduce  the  following: 
Definition  (Ref.  19):  The  system  described  by  Eqs.  (44)  Is  finite-time 
stable  (FTS)  In  the  interval  T  *  [tg,tf)  with  respect  to  (Y.8.t0,ll  l) 

for  every  trajectory  x^(t)  If  lXg(t0)l  <  y  Implies  BXg(t)l  <  8  for  every 

teT,  where  tt  i  denotes  the  Euclidean  norm. 

Moreover,  we  Introduce  the  notation 


B(a)  *  {xc;  ixci  <  a>,  V4-(t)  =  min  V(xc,t) 

(A1.A2) 

lxc»  -  a 

V4(t)  *  max  V(xc,t),  V®(t)  *  sup  V(xc,t) 

(A3.A4) 

"5c1  *  4  "5c"  =  a 


This  permits  us  to  state  the  following: 

Theorem  (Ref.  20):  The  system  described  by  Eqs.  (44)  Is  FTS  with 
respect  to  (Y,s,t0,tf ,11  II).  ris,  If  there  exists  a  function  V(xc,t), 
as  well  as  a  function  p(t)  integrable  on  T,  such  that: 

a.  V(xc,t)  <  p(t)  for  every  xceB(a)  and  te(t0,tf) 

b.  {  '  »(.)dT  <  vj(t)  -  »»(’>(to). 

We  choose 

V(xc,t)  *  xJPxc  (A5) 

Then,  taking  the  time  derivative  and  using  Eqs.  (44),  (46)  and  (47),  we 
obtain 

'J(Xc.t)  =  -  xJ(Qj  +  PGS_1GTP)xc  (A6) 

But, 

Q1  +  PGS_1GTP  >0,  teT  (A7) 

Hence,  Condition  a  of  the  theorem  Is  satisfied  by  choosing  0(t)  =  0,  so 
that  Condition  b  reduces  to 

vjW(t0)  S  vj(t),  teT  (A8) 
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Considering  Eq.  (A5),  Eq.  (A8)  yields 

xM[P(t0)lY2  *  max  xJ(t0)P(tQ)xc(t0)  <  min  xjp(t)xc  =  xm[P(t)]s2 

|XC»  *  Y  BXc8  »  a  (A9 

where  xM  and  xm  are  the  maximum  and  minimum  eigenvalues  of  the  argument 

matrix.  We  note  here  that  all  the  eigenvalues  of  P  are  real  and 
positive  In  T,  because  P  Is  symmetric  and  positive  definite. 


Introducing  llxc!  =  y  In  Eq.  (A9),  we  conclude  that 
|xc(t) If  <  8  (All) 


8  <  »Xc(t0)( 

Recalling  Eq.  (43a),  letting  tQ  *  0  and  introducing  the  notation 

XJIP1  *  min  xm(P(t)J,  teT 

we  obtain  _ _ 

A,[P(0,a)] 

ixr(t)l  <  lxr(0)i  /  — - e"  ,  teT 

-C  -C  /  XJIP(.)1 


/x„IP(t0)  1 
/ 


(A10) 

(A12) 

(A13) 


Inequality  (A13)  enables  us  to  state  that: 

I.  lx^(t ) B  approaches  zero  not  slower  than  e~at  and 

II.  In  most  cases,  the  ratio  of  the  eigenvalues  in  (A13)  does  not 
Increase  exponentially  with  a.  Hence,  Increasing  a  causes  x-(t)  to 
decay  to  zero  within  the  time  Interval  tf. 


Figure  1  -  The  Flexible  Spacecraft 
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Figure  2  -  Rigid-Body  Maneuver  (20°) 
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Figure  3  -  Output  Moments  From  Actuators  on  the  Beam  (20°) 
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Figure  4  -  Slope  Error  at  the  Tip  of  the  Beam  (180°) 
Four  Admissible  Functions  in  the  Model 


Figure  5  -  Slope  Error  at  the  Tip  of  the  Beam  (180°) 
Seven  Admissible  Functions  in  the  Model 
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Figure  6  -  Slope  Error  at  the  Tip  of  the  Beam  (20°) 
Four  Admissible  Functions  in  the  Model 
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Figure  7  -  Slope  Error  at  the  Tip  of  the  Beam  (20°) 
Seven  Admissible  Functions  in  the  Model 
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Abstract 

Undamped  distributed  structures  represent 
self-adjoint  systems,  admitting  real  eigenvalues 
and  real  orthogonal  eigenfunctions.  Control  of 
self-adjoint  systems  can  be  carried  out  conve¬ 
niently  by  modal  control.  Distributed  structures 
with  proportional  damping  possess  the  same  eigen¬ 
functions  as  the  corresponding  undamped  struc¬ 
tures,  so  that  the  same  modal  approach  can  be 
used  In  this  case  as  well.  Nonproportlonal 
damping  tends  to  destroy  the  self-adjointness  of 
the  system,  so  that  modal  control  Is  not  as  con¬ 
venient  as  for  undamped  structures.  If  damping  Is 
relatively  small,  however.  It  Is  possible  to  base 
the  control  design  on  the  self-adjoint  system  and 
still  obtain  satisfactory  control  performance. 


Introduction 

Structures  are  distributed-parameter  systems 
whose  motion  can  be  described  by  partial  differ¬ 
ential  equations  characterized  by  self-adjoint 

stiffness  operators1.  In  the  absence  of  damping, 
this  Implies  that  the  structures  themselves  are 
self-adjoint,  which  Implies  further  that  the 
eigenfunctions  are  orthogonal.  In  general, 
damping  tends  to  destroy  the  sel f-ad jointness 
property.  One  notable  exception  is  proportional 
damping,  for  which  the  eigenfunctions  are  the 
same  as  for  the  corresponding  undamped  struc¬ 
tures,  so  that  the  orthogonality  is  preserved. 

The  self-adjointness  property  of  distributed 
structures  is  very  important  in  control  design. 
Indeed,  for  undamped  structures  it  Is  possible  to 
take  advantage  of  the  eigenfunctions  orthogonal¬ 
ity  and  transform  the  partial  differential  equa¬ 
tion  Into  a  set  of  Independent  second-order  ordi¬ 
nary  equations,  known  as  modal  equations.  Then, 
controls  can  be  designed  for  each  modal  equation 
independently ,  making  It  possible  to  target  Indi¬ 
vidual  modes  for  control2-5.  The  same  is  true 
for  systems  with  proportional  damping.  In  the 
general  case  of  feedback  control  of  non-sel f- 
adjolnt  distributed  structures,  full  decoupling 
of  the  distributed  system  may  not  be  possible,  so 
that  modal  control  loses  some  of  Its  appeal. 

In  many  cases,  the  factors  rendering  the 
structures  non-sel f-adjolnt  tend  to  be  suf¬ 
ficiently  small  that  they  can  be  regarded  as 
perturbations  on  the  self-adjoint  system.  Then, 


the  question  arises  as  to  the  possibility  of 
basing  the  controls  on  the  self-adjoint  system. 
This  Is  the  approach  explored  in  this  paper.  A 
particular  question  to  be  answered  is  how  large 
the  factors  rendering  the  system  non-self-adjoint 
can  be  without  invalidating  the  approach.  This 
question  Is  intimately  related  to  that  of  robust¬ 
ness,  i.e.,  how  sensitive  Is  the  control  design 
to  changes  in  the  system  parameters,  where  these 
changes  are  the  ones  responsible  for  the  non- 
sel  f-adjointness  of  the  system.  One  way  of  look¬ 
ing  at  the  problem  Is  by  examining  the  closed- 
loop  poles.  In  particular,  the  question  Is 
whether  the  factors  rendering  the  structure  non- 
self-adjolnt  can  push  the  closed-loop  poles  far 
away  from  the  closed-loop  poles  obtained  from  the 
control  design  based  on  the  self-adjoint  struc¬ 
ture,  thus  Impairing  the  control  system  perfor¬ 
mance. 

Quite  often,  control  of  structures  Is 
carried  out  as  If  the  structures  were  undamped. 
Inherent  damping  in  the  structure,  however, 
causes  the  actual  closed-loop  poles  to  differ 
from  the  modeled  closed -loop  poles.  It  is  shown 
In  this  paper  that  if  the  control  system  design 
takes  advantage  of  the  self-adjointness  proper¬ 
ties  of  structures,  then  the  sensitivity  of  the 
closed-loop  poles  to  the  addition  of  small 
general  viscous  damping  can  be  reduced  to  a 
simple  algebraic  expression.  In  fact,  the 
closed-loop  eigenvalue  sensitivity  turns  out  to 
be  a  function  of  the  diagonal  damping  coef¬ 
ficients,  so  that  the  eigenvalues  are  insensitive 
to  damping  that  might  cause  the  system  to  become 
non-self-adjoint.  A  control  system  design  which 
takes  advantage  of  the  self-adjoint  property  In 
structures  is  Independent  modal -space  control 

(IMSC)2-5. 

This  paper  examines  the  circumstances  under 
which  controls  based  on  a  self-adjoint  structure 
yield  satisfactory  performance  when  applied  to 
the  corresponding  non-self-adjoint  structure.  To 
this  end,  the  partial  differential  equation  of 
motion  for  a  distributed  structure  Is  first  In¬ 
troduced.  Then,  a  perturbation  technique  is  used 
to  compute  the  sensitivity  of  the  closed-loop 
poles  to  the  addition  of  small  viscous  damping  in 
the  system.  The  paper  follows  with  an  investi¬ 
gation  of  the  sensitivity  of  the  closed-loop 
poles  when  advantage  is  taken  of  the  self-adjoint 
properties  to  design  independent  modal -space 
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control  for  the  structure.  A  numerical  example 
Is  presented  In  which  the  control  system  perfor¬ 
mance  Is  examined  using  varying  amounts  of 
damping.  The  numerical  example  indicates  in  a 
heuristic  manner  how  large  the  nonproportional 
damping  can  be  without  degrading  control  system 
performance. 

It  Is  perhaps  appropriate  to  mention  another 
work  using  a  perturbation  approach  to  the  control 

of  structures6.  The  structure  considered  in  Ref. 
6  1$  undamped  and  the  control  Is  carried  out  by  a 
low-authority  controller,  i.e.,  a  controller  pro¬ 
viding  sufficiently  small  damping  that  the 
modeled  closed-loop  poles  could  be  computed  from 
the  open-loop  elgensolutlon  on  the  basis  of  a 

perturbation  technique*.  In  contrast,  the 
structure  considered  here  Is  inherently  damped, 
but  the  damping  Is  sufficiently  small  that  the 
controller  can  be  designed  by  ignoring  the 
damping.  Then,  the  perturbation  approach  is  used 
to  calculate  the  shift  In  the  undamped  system 
closed-loop  poles  caused  by  the  Inherent  damping. 


Equations  of  Motion 

The  motion  of  a  distributed  structure  can  be 

derived  using  the  extended  Hamilton's  principle* 
and  is  given  by  the  partial  differential  equation 

u(P,t)  ♦  u(P,t)  +  m(P)u(P,t)  *  f(P,t),  PeD 

(1) 

where  u(P,t)  is  the  displacement  at  point  P  In 
the  domain  D  of  the  structure,  is  a  self- 
adjoint  differential  operator  of  order  2p  expres¬ 
sing  the  system  stiffness.  Is  a  viscous  damping 
operator,  m(P)  is  a  scalar  function  representing 
the  mass  distribution  and  f(P,t)  Is  the  external 
force  density.  The  displacement  u(P,t)  Is 
subject  to  the  boundary  conditions 

Biu(P,t)  *  0,  1  *  1,2,. ...p;  PeS  (2) 

where  B,-  are  differential  operators  of  maximum 
order  2p-l  and  S  defines  the  boundary  of  0. 

We  consider  the  eigenvalue  problem  associ¬ 
ated  with  the  undamped  self-adjoint  structure 
gi  ven  by 

♦r(p)  •  V"(P)V(P)-  r  *  PeD  <3> 

and  the  boundary  conditions 

♦  (P)  ■  0,  1  ■  1,2,. ...p;  r  *  1,2,...;  PeS 

(4) 

The  solution  to  Eqs.  (3)  and  (4)  consists  of  a 
denumerably  Infinite  set  of  real  nonnegative 
eigenvalues  *r  and  associated  real  eigenfunctions 

♦r(P).  The  eigenvalues  are  related  to  the 

natural  frequencies  of  undamped  vibration 
2 

by  *r  ■  «r  (r  *  1,2,...).  For  convenience,  we 
order  the  eigenvalues  so  that  Xj  <  x?  <  ...  . 


The  eigenfunctions  can  be  normalized  so  as  to 
satisfy 

/0  mer(p)as(p)  dO.  /D  ♦P(P)  ♦S(P)  do  «  uJsrs, 

r,s  *  1,2, ...  (5) 

where  4ps  is  the  Kronecker  delta  function. 

Equations  (5)  represent  the  orthonormality 
conditions  for  self-adjoint  systems. 

We  propose  to  express  the  solution  of  Eqs. 
(1)  and  (2)  In  terms  of  the  Infinite  series 

u(P.t)  -  I  ♦r(P)up(t)  (6) 

r*l  r  r 

where  up(t)  are  time-dependent  generalized 

coordinates.  Introducing  Eq.  (6)  into  (1), 
multiplying  the  result  by  integrating  over 

the  domain  0  and  considering  Eqs.  (5),  we  obtain 
ur(t)  +  l  crsus(t)  +  u>pur(t)  ■  fp(t), 

r  ■  1.2,...  (7) 

where 

Crs’JoV  *s  dD*  r,s  *  1.2,...  (8) 

are  viscous  damping  coefficients  and 

fr(t)  -  /„  ♦r(P)f (P,t)  dO,  r  •  1.2,...  (9) 

are  generalized  control  forces. 

Equations  (7)  represent  an  infinite  set  of 
ordinary  differential  equations  expressing  the 
motion  of  a  distributed  structure.  In  the 
absence  of  damping,  crs  »  0  (r.s  *  1,2,...),  the 

system  Is  self-adjoint  and  Eqs.  (7)  represent  a 
decoupled  set.  However,  damping  tends  to  prevent 
the  decoupling  and  hence  to  destroy  the  self¬ 
adjointness  property.  In  the  case  of  non-self- 
adjoint  systems  It  is  necessary  to  first  define 

and  then  solve  an  adjoint  eigenvalue  problem*. 
Only  in  special  cases  does  the  operator  have 
♦r  (r  *  1,2,...)  as  its  eigenfunctions.  In 

these  cases,  crs  *  0  (r*s).  Note  that  when  the 
operator  admits  ar  (r  «  1,2,...)  as  its  eigen¬ 
functions,  the  normal  modes  of  the  undamped 
structure  are  Identical  to  the  normal  modes  of 

the  damped  structure*. 

Equations  (7)  can  be  expressed  in  the  vector 

form 


u(t)  +  Cu(t)  +  Au(t)  *  f(t) 

(10) 

where 

u(t)  »  Cuj(t)  u2(t)  ...] 

(11) 

Is  an  Infinite-dimensional  vector 
coordinates. 

of  generalized 

f(t)  -  Cfj(t)  f2(t)  ...] 

(12) 

Is  an  Infinite-dimensional  vector  of  generalized 
forces , 

A  *  diag  [u2  u2  ...]  (13) 

is  the  ««•  matrix  of  natural  frequencies  of  the 
undamped  structure  squared  and  C  »  [c„]  is  the 
damping  matrix. 

In  general,  the  damping  matrix  C  is  fully 
populated.  In  this  case,  to  obtain  a  solution, 
it  Is  necessary  to  cast  Eq.  (10)  in  the  state- 
space.  To  this  end,  we  introduce  the  2«- 

dimensional  state  vector  x(t)  *  [uT(t)  uT(t)]T. 

Then,  Eq.  (10)  can  be  written  in  the  state  vector 
form 

x  *  A*x  +  8*f  (14) 


-L-il .  b*  *  r?" 

!  -  cj  |_i  _ 


are  coefficient  matrices,  in  which  I  is  an  «x. 
identity  matrix. 


Sensitivy  Analysis 


We  consider  the  case  In  which  damping  Is 
sufficiently  small  that  the  control  system  design 
can  be  based  on  the  undamped  model ,  where  the 
latter  is  obtained  by  letting  C  *  0  in  Eq.  (10). 
In  general,  the  feedback  control  recouples  the 
undamped  equations.  Indeed,  the  control  law  is 
given  by 

f(t)  «  -  Gu  -  Hu  (16) 

where  in  general  G,  H  are  fully  populated  «•- 
dimensional  control  gain  matrices  with  entries 
gPS,  hrs  (r,s  =  1,2,...),  respectively.  The  net 

effect  of  the  feedback  control  is  to  recouple  the 
undamped  equations  of  motion,  Eq.  (10),  through 
the  nonzero  off-diagonal  terms  in  the  gain 
matrices  G  and  H.  Substituting  Eq.  (16)  into  Eq. 
(10)  with  C  *  0,  we  obtain  the  closed-loop  state 
equations  of  motion  in  the  vector  form 

x  «  AgX  (17) 


The  eigenvalue  problems  associated  with  Eq.  (17) 
are 

A0 ~0r  *  X0r~0r’  A0^0s  *  X0s~0s*  r,S  *  1’2”'* 


(19a, b) 


where  x0r  (r  -  1.?....)  are  the  closed-loop  poles 

of  the  undamped  system.  We  note  that,  because  Aq 

is  not  symmetric,  we  must  consider  both  the  right 
and  left  eigenvectors  of  Aq,  denoted  by  u0f>  and 


Vqs,  respectively.  The  two  sets  of  eigenvectors 
~0r*  ~0$  are  t>1ortho9onal  and  can  be  normalized 
so  as  to  satisfy  the  biorthonormal ity  conditions* 


KOs^Or 


5rs*  ~OsAO^Or  *  x0r6rs* 


r,s  *  1,2,...  (20a, b) 


In  general,  the  eigenvalues  AQr  (r  »  1,2,...)  and 

eigenvectors  un  ,  vn  (r,s  *  1,2,...)  are  complex 
quantities.  r  ~us 

By  virtue  of  the  small -damping  assumption, 
the  open-loop  actual  distributed  structure  is 
nearly  self-adjoint,  and  differs  from  the  open- 
loop  undamped  structure  by  a  small  perturbation. 
Of  course,  the  perturbation  arises  from  the  fact 
that  the  actual  distributed  structure  contains 
small  viscous  damping.  The  actual  closed-loop 
equations  of  motion  are  given  by 

x  =  Ax  (21) 


.  i  ,1 

-  A  -  G  |  -  C  -  H 


The  eigenvalue  problems  associated  with  Eq.  (21) 
have  the  same  form  as  Eqs.  (19),  except  that  the 
subscript  0  is  not  present.  Moreover,  the  two 
sets  of  eigenvectors  up,  vr  can  be  normalized  so 

as  to  satisfy  the  biothonormal ity  conditions 
given  by  Eqs.  (20)  with  the  subscript  0  omitted. 

When  damping  in  the  actual  distributed 
structure  is  small,  i.e.,  when  crs  =  0(c)  (r,s  * 

1,2,...),  the  matrix  A  can  be  expressed  in  the 
perturbed  form 

A  =  Aq  +  Az  (23) 

where  Aq  is  given  by  Eq.  (18)  and  Aj  is  a  first- 
order  perturbation  matrix  given  by 


o  _ o 

0  j-  C 


Because  A  is  obtained  from  A0  by  a  small  pertur¬ 
bation,  we  assume  that  the  eigensol utions  can  be 
written  in  the  form 


V  ’  X0r  +  Xlr’ 


iiOr  +  ~Ir* 


*0r  *  ~lr* 


1 , 2 ,  •  • , 


r  «  1,2,... 


1,2,... 


It  can  be  shown  that  the  eigenvalue  and  eigen¬ 
vector  perturbations  have  the  expressions* 


~0rAl~0r*  r 


1,2.... 
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(26b) 
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(26c) 


Equations  (26a)  can  be  regarded  as  providing  the 
sensitivity  of  the  closed-loop  poles  to  small 
viscous  damping  in  the  distributed  structure. 
Indeed,  dividing  both  sides  of  Eqs.  (26a)  by  c1 j , 

we  obtain  the  first-order  eigenvalue  sensitivity 
coefficient  3Ar/3c..j,  which  is  a  measure  of  the 

change  in  the  closed-loop  eigenvalue  Af  due  to 

the  damping  coefficient  c 


i  j  ’ 

closed-loop  eigenvectors  uQr  and  vQr  (r 


In  general ,  the 

1. 

2,...)  are  fully  populated  so  that  the  perturba¬ 
tions  Alf.,  ulr>  vlr  (r  »  1,2,...)  are  sensitive 

to  all  the  damping  terms  cps  (r,s  »  1,2,...). 


Sensitivity  Using  Natural  Control 

Natural  control  is  characterized  by  the 
preservation  of  the  natural  coordinates  in  the 
closed-loop  system,  which  can  be  traced  to  the 
fact  that  the  open-loop  eigenfunctions  are  iden¬ 
tical  to  the  closed-loop  eigenfunctions.  It  can 
be  shown  that  natural  control  is  the  globally 
optimal  solution  to  the  control  problem  for  self- 

adjoint  distributed  parameter  systems3.  Natural 
control  is  obtained  using  the  Independent  modal 

space  control  (IMSC)  method^.  The  IMSC  control 
gain  matrices  G  and  H  are  diagonal  with  entries 
equal  to  gr  (r  »  1,2,...)  and  hr  (r  «  1,2,...), 

respectively.  Hence,  natural  control  can  be  re¬ 
garded  as  providing  proportional  viscous  damping 
and  proportional  stiffness  augmentation. 

Recognizing  that  the  right  eigenvectors  have 
the  form  uQr  -  [aJr '  V^]1,  *here  30r  <r  "  l> 
2,...)  are  —dimensional  configuration  eigen¬ 
vectors,  Eq.  (19a)  can  be  written  in  the  explicit 
form 


0 

I 

fteJl 

i- 

St 

1 _ 

-  A  -  G 

-  H 

Ji 

O 

o 

*< 

_ 1 

*  *0r 

X0rSdr 

(27) 


It  can  be  shown  that  when  G  and  H  are  diagonal 
matrices,  (jQr  *  er,  where  ep  is  a  standard  unit 

vector  (with  all  its  components  equal  to  zero 
with  the  exception  of  the  rth  component  which  is 
equal  to  one).  Then,  from  Eqs.  (19b)  and  (20a), 
we  obtain 


”  Sr" 

V  -  -L 

(ur  +  gr] Sr 

HOr  * 

x0rSr 

Lm  — 

*  ~0r  ap 

_  *  x0rSr  _ 

1.2.. 


(28a, b) 


Hence,  in  the  case  of  IMSC,  characterized  by 
diagonal  control  gain  matrices  G  and  H  in  Eq. 
(16),  each  of  the  closed-loop  eigenvectors  u0r, 

Vqp  (r  ■  1,2,...)  associated  with  the  undamped 

structure  contains  only  two  nonzero  terms.  The 
closed-loop  eigenvalues  A0r  and  the  normal izatlon 

constants  ar  (r  «  1,2,...)  have  the  expressions 
x0r  *  *  7  {hr  +  1[4(u?  +  9r>  *  h?]1/2}  (29a) 

°r  *  “p  +  9r  -  x0r  (29b) 

Equations  (28)  and  (29)  specify  only  one  half  of 
the  eigenvalues  and  left  and  right  eigenvectors. 
The  complex  conjugates  comprise  the  other  half. 


Inserting  Eqs.  (28)  into  (26a),  we  can  write 


2 

where  the  real  part  of  the  quantity  AQr/ar  (r  * 

1,2,...)  is  always  negative.  The  perturbation  in 
the  right  eigenvectors  is  obtained  from  Eq.  (26b) 
and  is  given  by 


1  *  4kr  X0r  x0k 


Ur 
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Or 


'Ok 


'  ckr  -Ok  * 

P  *  1,2,...  (31) 


Examining  Eq.  (30),  for  a  given  closed-loop 
system.  It  is  apparent  that  the  perturbations  in 
the  closed-loop  eigenvalues  vary  only  as  a 
function  of  the  diagonal  terms  in  the  matrix  C. 
Hence,  damping  perturbations  in  which  crp  >  0 

(r  *  1,2,...)  cannot  shift  the  actual  structure 
closed-loop  poles  to  the  right  in  the  complex 
plane  relative  to  the  corresponding  closed-loop 
poles  of  the  undamped  structure.  Note  that  this 
conclusion  was  reached  on  the  basis  of  natural 
control.  The  off-diagonal  terms  in  the  damping 
matrix  are  the  terms  inducing  the  actual  dis¬ 
tributed  structure  to  vibrate  in  a  linear  combi¬ 
nation  of  complex  modes,  which  are  different  from 
the  undamped  normal  modes.  However,  when  the 
damping  terms  are  small,  as  is  the  case  here,  the 
closed-loop  poles  are  sensitive  only  to  changes 
in  the  damping  matrix  which  do  not  Induce  non- 
self-adjoint  properties.  Moreover,  when  the 
damping  matrix  C  is  diagonally-dominant,  the 
perturbations  in  the  eigenvectors  given  by  Eqs. 
(31)  are  of  second-order  in  magnitude,  and  the 
closed-loop  system  tends  to  vibrate  essentially 
as  a  linear  combination  of  its  normal  modes  only. 


Numerical  Example 

To  illustrate  the  sensitivity  of  the  closed- 
loop  poles  and  control  system  performance  to 
small  viscous  damping,  we  consider  the  control  of 


the  bending  vibration  of  a  simply-supported 
beam.  Assuming  that  the  beam  has  uniform  mass 

and  stiffness  properties,  Eq.  (1)  becomes7 


Note  that  with  6  »  0,  we  have  proportional 
damping,  in  which  case  we  obtain  c„ 
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r  rs 


a4u  s2 

El  ^-jf  +  ^-r  [cl (x) 

3X  3xc 


S3u  i 
~T~. 
3x  at 


♦  N  *  f(x,t), 
3t‘ 

0  <  x  <  L  (32) 


where  we  note  that  the  damping  operator  in  Eq. 

2  2 

(32),  given  by  «  [cl(x)  ^~w],  provides 
3x  3x 

strain  rate  damping,  which  is  a  form  of  viscous 
damping.  The  displacement  u  is  subject  to  the 
boundary  conditions 


u(0,t) 


a2u(o,t) 


u(L,t) 


?lu-lL-ali  ,  o  (33) 

3x“" 


For  simplicity,  we  choose  M  •  1,  El  »  1  and  L  *  5. 

It  can  be  shown  that  the  solution  of  the 
associated  eigenvalue  problem,  Eq.  (3),  for  the 
undamped  case,  i.e.,  for  c I ( x)  *  0,  consists  of 

the  eigenfunctions  and  natural  frequencies®, 

♦r(x)  -  (£)1/2s1n  (JfL).  r  «  1,2,...  (34a) 


wr  *  (p-)2.  r  «  1,2,...  (34b) 

respectively.  For  the  damped  case,  we  choose  the 
strain  rate  damping  distribution  to  be  given  by 

cl(x)  *  a  *  8x(L  -  x),  0  <  x  <  L  (35) 


where  a  and  a  are  constants.  Note  that  when  8  * 
0,  the  strain  rate  damping  is  uniform,  which  re¬ 
sults  in  the  case  of  proportional  damping,  as  in 
this  case  the  damping  distribution  Is  propor¬ 
tional  to  the  stiffness  distribution.  However, 
the  presence  of  a  nonzero  8  destroys  the  self¬ 
adjointness  property,  producing  nonzero  off- 
diagonal  entries  in  the  damping  matrix  C.  From 
Eq.  (8),  we  obtain  upon  integrations  by  parts 

2  <j2. 

Crc  *  In  ♦r(x^  (£•*  +  6x(L  *  —r-)  <*x 

rs  u  dx‘  dx‘ 


,  d24.(x)  d2*  (x) 

/o  [a  +  8x(L  -  x)]  — - -j —  dx 


dx 


dx 


(36) 


With  the  eigenfunctions  given  by  Eq.  (34a),  the 
damping  coefficients  In  Eq.  (36)  can  be  computed 
in  closed-form  and  are  given  by 
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..2.2  2  ,  . ,r-s  ,  ,  . >r+s  , 
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(37) 


r, s  »  1,2, 


We  consider  controlling  a  subset  of  the 
modes  using  IMSC,  as  controlling  the  entire 
infinity  of  modes  would  require  a  distributed 

actuator2.  More  specifically,  the  first  three 
modes  are  to  be  controlled  using  3  discrete 
actuators  at  locations  x^  *  1L/4  (1  «  1,2,3) 

together  with  modal  filters2,  which  are  used  for 
the  extraction  of  the  modal  coordinates.  To 
determine  the  control  gain  matrices  G  and  H  In 
Eq.  (16),  the  following  performance  Index  is 
minimized 

J  *  J  /q  +  V?(t)  +  VrU)3  «  (38) 

Note  that  the  first  two  terms  in  the  integrals  in 
/r  Eq.  (38)  represent  modal  kinetic  and  potential 
energies,  respect! vely,  and  the  third  term  repre¬ 
sents  the  modal  control  effort,  all  for  the  first 
three  modes.  The  minimization  yields  the  modal 

gains2 

9r  =  -  a2  +  <ur(uZ  ♦  1/Rr)1/2,  r  *  1,2,3  (39a) 

hr  *  C-  *  1/Rr  ♦  aap(wj  +  1/Rr)l/2]1/Z, 

r  -  1,2,3  (39b) 

The  control  law  given  by  Eq.  (16)  with  the  con¬ 
trol  gains  provided  by  Eqs.  (39)  represents 
linear  optimal  IMSC. 

To  show  the  effects  of  damping,  we  investi¬ 
gate  the  3  cases,  a  «  B  *  0,  a  =■  b  =  0.0001  and 
a  *  8  *  0.001,  denoted  by  I,  II  and  III,  respec¬ 
tively.  Table  1  displays  the  damping  matrix  with 
entries  given  by  Eq.  (37)  with  a  *  8  *  1,  where 
the  actual  distributed  structure  is  simulated  by 
truncating  the  infinite-dimensional  system  by 
retaining  8  terms  in  series  (6).  Table  2  pro¬ 
vides  the  closed-loop  poles  for  the  3  cases  with 
Rr  *  1.0  (r  *  1,2,3).  Examining  Figure  1,  we 

note  that  damping  shifts  the  closed-loop  poles  to 
the  left  in  the  complex  plane  as  the  level  of 
damping  increases.  Morever,  the  shift  in  the 
closed-loop  poles  for  the  controlled  modes,  i.e., 
the  lowest  3  modes,  can  be  computed  using  Eqs. 
(29)  and  (30),  where  the  natural  frequencies  u 
(r  *  1,2,...)  are  given  by  Eq.  (34b)  and  the 
control  gains  gr,  hr  (r  «  1,2,3)  are  computed 

using  Eqs.  (39)  with  Rr  =  1.0  (r  *  1,2,3). 

To  examine  the  control  system  performance, 
we  considered  Cases  I,  II  and  III,  as  above, 
except  that  we  increased  the  control  effort  by 
decreasing  Rr  in  Eq.  (38).  To  this  end,  we  used 

Rr  ■  0.01  (r  •  1,2,3),  which  has  the  effect  of 

increasing  the  gains,  as  can  be  seen  from  Eqs. 
(39).  Figures  2  and  3  plot  the  performance  index 
J(t)  given  by  Eq.  (38)  for  the  values  Rr  »  1.0 

and  Rr  »  0.01  (r  »  1,2,3),  respectively,  for  the 

above  three  cases.  It  is  evident  from  Figure  2 
that  for  Rr  *  1.0  the  effect  of  damping  is  to 


;r  *  S 


enhance  the  control  system  performance,  as  the 
performance  Index  J(t)  has  the  values  893.1, 

883.7  and  813.1  at  t  «  5  s  corresponding  to  the 
Cases  I,  II  and  III,  respectively.  Figure  3, 
however,  displays  the  opposite  effect,  with  the 
performance  index  J(t)  for  Rr  *  0.01  having  the 

values  250.3,  250.4  and  252.3  at  t  *  5  s  corres¬ 
ponding  to  the  Cases  I,  II  and  III,  respectively. 
Hence,  damping  may  enhance  or  degrade  control 
system  performance  depending  on  the  value  of  the 
gains.  Indeed,  as  the  gains  Increase,  larger 
amounts  of  energy  are  transferred  back  and  forth 
between  the  controlled  and  uncontrolled  modes, 
which  are  coupled  due  to  the  presence  of  damping. 

This  effect  is  known  as  control  spillover®.  How¬ 
ever,  we  also  conclude  that  the  damping  does  not 
alter  control  system  performance  significantly. 


Conclusions 

Quite  often,  control  of  structures  is 
carried  out  as  if  the  structures  were  undamped. 
Undamped  structures  are  self-adjoint.  Inherent 
structural  damping  tends  to  destroy  the  self- 
adjointness  properties  and  can  degrade  control 
system  performance.  However,  designing  controls 
on  the  basis  of  an  undamped  structures  has  many 
advantages.  One  of  these  advantages  is  that  is 
permits  the  use  of  IMSC.  It  turns  out  that  in 
this  case  the  perturbation  in  the  modeled  closed- 
loop  poles  due  to  small  damping  does  not  cause  a 
dramatic  shift  in  the  closed-loop  poles  nor  does 
it  change  the  control  system  performance  very 
much.  The  sensitivity  study  presented  here  indi¬ 
cates  that  damping  shifts  the  actual  closed-loop 
poles  to  the  left  of  the  modeled  closed-loop 
poles  in  the  complex  plane.  Moreover,  the 
presence  of  small  damping  in  the  actual  dis¬ 
tributed  structure  does  not  affect  the  control 
system  performance  significantly  in  the  case  of 
IMSC. 
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